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Abstract 

When modelling metapopulation dynamics, the influence of a sin- 
gle patch on the metapopulation depends on the number of individu- 
als in the patch. Since the population size has no natural upper limit, 
this leads to systems in which there are countably infinitely many 
possible types of individual. Analogous considerations apply in the 
transmission of parasitic diseases. In this paper, we prove a law of 
large numbers for quite general systems of this kind, together with 
a rather sharp bound on the rate of convergence in an appropriately 
chosen weighted £i norm. 
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1 Introduction 

There are many biological systems that consist of entities that differ in their 
influence according to the number of active elements associated with them. 
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and can be divided into types accordingly. In parasitic diseases (Barbour & 
Kafctzaki 1993, Luchsinger 2001a,b, Kretzschmar 1993), the infectivity of a 
host depends on the number of parasites that it carries; in metapopulations, 
the migration pressure exerted by a patch is related to the number of its 
inhabitants (Arrigoni 2003); the behaviour of a cell may depend on the num- 
ber of copies of a particular gene that it contains (Kimmel & Axelrod 2002, 
Chapter 7); and so on. In none of these examples is there a natural upper 
limit to the number of associated elements, so that the natural setting for 
a mathematical model is one in which there are countably infinitely many 
possible types of individual. In addition, transition rates typically increase 
with the number of associated elements in the system — for instance, each 
parasite has an individual death rate, so that the overall death rate of par- 
asites grows at least as fast as the number of parasites — and this leads 
to processes with unbounded transition rates. This paper is concerned with 
approximations to density dependent Markov models of this kind, when the 
typical population size N becomes large. 

In density dependent Markov population processes with only finitely 
many types of individual, a law of large numbers approximation, in the form 
of a system of ordinary differential equations, was established by Kurtz (1970), 
together with a diffusion approximation (Kurtz, 1971). In the infinite di- 
mensional case, the law of large numbers was proved for some specific mod- 
els (Barbour & Kafetzaki 1993, Luchsinger 2001b, Arrigoni 2003, see also 
Leonard 1990), using individually tailored methods. A more general result 
was then given by Eibeck & Wagner (2003). In Barbour & Luczak (2008), 
the law of large numbers was strengthened by the addition of an error bound 
in ii that is close to optimal order in N. Their argument makes use of an 
intermediate approximation involving an independent particles process, for 
which the law of large numbers is relatively easy to analyse. This process is 
then shown to be sufficiently close to the interacting process of actual inter- 
est, by means of a coupling argument. However, the generality of the results 
obtained is limited by the simple structure of the intermediate process, and 
the model of Arrigoni (2003), for instance, lies outside their scope. 

In this paper, we develop an entirely different approach, which circum- 
vents the need for an intermediate approximation, enabling a much wider 
class of models to be addressed. The setting is that of famihes of Markov 
population processes := (Xiv(i), t > 0), N > 1, taking values in the 
countable space A'+ := {X e Z^+; ^^>qX'^ < oo}. Each component repre- 
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sents the number of individuals of a particular type, and there are countably 
many types possible; however, at any given time, there are only finitely 
many individuals in the system. The process evolves as a Markov process 
with state-dependent transitions 

X -> X + J at rate NajiN'^X), X e X+, J e J, (1.1) 

where each jump is of bounded influence, in the sense that 

J7 C {X G Z^+; ^ |X™| < J, < cx)}, for some fixed J, < oo, (1.2) 

m>0 

SO that the number of individuals affected is uniformly bounded. Density 
dependence is reflected in the fact that the arguments of the functions aj 
are counts normalised by the 'typical size' X. Writing 7^ := -R_,_^, the func- 
tions af. IZ — )► i?+ are assumed to satisfy 

5^aj(0 < oo, ee7^o, (1.3) 

where TZq := G 7^: = for all but finitely many i}; this assumption 
implies that the processes X^v are pure jump processes, at least for some 
non-zero length of time. To prevent the paths leaving Xj^, we also assume 
that Ji > —1 for each /, and that a;j(^) = if = for any J E J such 
that J' = —1. Some remarks on the consequences of allowing transitions J 
with J' < —2 for some / are made at the end of Section |H 

The law of large numbers is then formally expressed in terms of the system 
of deterministic equations 

I = 5^Ja,(0 =: i^o(0, (1-4) 

to be understood componentwise for those ^ e71 such that 
^|JVj(0 < oo> for all />0, 

thus by assumption including 71q. Here, the quantity Fq represents the in- 
finitesimal average drift of the components of the random process. However, 
in this generality, it is not even immediately clear that equations (11. 4p have 
a solution. 
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In order to make progress, it is assumed that the unbounded components 
in the transition rates can be assimilated into a hnear part, in the sense 
that Fq can be written in the form 

Fo(0 = Ae + F(0, (1.5) 

again to be understood componentwise, where A is a constant Z+ x Z+ 
matrix. These equations are then treated as a perturbed hnear system 
(Pazy 1983, Chapter 6). Under suitable assumptions on A, there exists a 
measure fi on Z+, defining a weighted ii norm || ■ ||^ on TZ, and a strongly 
II • II ^-continuous semigroup {R{t), t > 0} of transition matrices having point- 
wise derivative -R'(O) = A. If F is locally || ■ ||^-Lipschitz and ||a;(0)||^ < oo, 
this suggests using the solution x of the integral equation 

x{t) = R{t)x{0)+ [ R{t~ s)F{x{s))ds (1.6) 
Jo 

as an approximation to xn '■= N~^Xn, instead of solving the deterministic 
equations (11. 4p directly. We go on to show that the solution Xn of the 
stochastic system can be expressed using a formula similar to (11. 6p . which 
has an additional stochastic component in the perturbation: 

XN{t) = R{t)xN{0)+ [ R{t~ s)F{xN{s))ds + mN{t), (1.7) 

where 







rt 

mN{t) := I R{t~ s)dmN{s), (1.8) 





and is the local martingale given by 



niNit) := XAr(t) - XAr(O) - / Fo{xNis))ds. (1.9) 

Jo 

The quantity m^r can be expected to be small, at least componentwise, under 
reasonable conditions. 

To obtain tight control over in all components simultaneously, suf- 
ficient to ensure that supo<s<i ||^Ar(s)||^ is small, we derive Chernoff-like 
bounds on the deviations of the most significant components, with the help of 
a family of exponential martingales. The remaining components are treated 
using some general a priori bounds on the behaviour of the stochastic system. 
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This allows us to take the difference between the stochastic and deterministic 
equations fll.7p and f ll.6p . after which a Gronwall argument can be carried 
through, leading to the desired approximation. 

The main result, Theorem \A.7\ guarantees an approximation error of or- 
der 0{N~^/'^^/[ogN) in the weighted ii metric || ■ ||^, except on an event of 
probability of order 0(A^~^ log A^). More precisely, for each T > 0, there 
exist constants K^\k!^\k!^^ such that, for large enough, if 

\\N-'Xr,{0)-xm,<K?^^, 

then 

p(^=up^||iV-'A-„(i) - .(i)IU > < (1.10) 

The error bound is sharper, by a factor of log N, than that given in Barbour & 
Luczak (2008), and the theorem is applicable to a much wider class of models. 
However, the method of proof involves moment arguments, which require 
somewhat stronger assumptions on the initial state of the system, and, in 
models such as that of Barbour & Kafetzaki (1993), on the choice of infection 
distributions allowed. The conditions under which the theorem holds can be 
divided into three categories: growth conditions on the transition rates, so 
that the a priori bounds, which have the character of moment bounds, can 
be established; conditions on the matrix A, sufficient to limit the growth of 
the semigroup R, and (together with the properties of F) to determine the 
weights defining the metric in which the approximation is to be carried out; 
and conditions on the initial state of the system. The a priori bounds are 
derived in Section |21 the semigroup analysis is conducted in Section [31 and 
the approximation proper is carried out in Section |H The paper concludes 
in Section [5] with some examples. 

The form (11.81) of the stochastic component rhiq{t) in (ll.7p is very simi- 
lar to that of a key element in the analysis of stochastic partial differential 
equations; see, for example. Chow (2007, Section 6.6). The SPDE arguments 
used for its control are however typically conducted in a Hilbert space con- 
text. Our setting is quite different in nature, and it does not seem clear how 
to translate the SPDE methods into our context. 
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2 A priori bounds 

We begin by imposing further conditions on the transition rates of the pro- 
cess Xn, sufficient to constrain its paths to bounded subsets of dur- 
ing finite time intervals, and in particular to ensure that only finitely many 
jumps can occur in finite time. The conditions that follow have the fiavour 
of moment conditions on the jump distributions. Since the index j G Z+ is 
symbolic in nature, we start by fixing an z/ G 7^, such that refiects in 
some sense the 'size' of j, with most indices being 'large': 

> 1 for all j > and lim z/(j) = oo. (2.1) 

We then define the analogues of higher empirical moments using the quanti- 
ties Uj. G 71, defined by z/j,(j) := i^{jY, r > 0, setting 

Sr{x) := ''^^h'r{j)x^ = X^Ur, X E TZq, (2.2) 

i>o 

where, for x G TZq and y E 71, x'^y := J2i>o^iyi- particular, for X G 
So{X) = Note that, because of f l2.ip . for any r > 1, 

#{X G X+: Sr{X) <K}<oo for all K > 0. (2.3) 

To formulate the conditions that limit the growth of the empirical moments 
of X^it) with t, we also define 

Urix) := ^aj{x)J^iyr; Vrix) := ^ aj(x)( J^z/^)^ x E 71. (2.4) 

The assumptions that we shall need are then as follows. 

Assumption 2.1 There exists a v satisfying f l2.ip and r^lx,r^lx > 1 such 
that, for all X G Xj^, 

J2MN-'X)\j''ur.\ < OO, 0<r<r«„ (2.5) 

the case r = following from (11. 2p and (11. 3p ." furthermore, for some non- 
negative constants kri, the inequalities 

Uo{x) < koiSo{x) + ko4, 

Ui{x) < knSi{x) + ku, (2.6) 

Ur{x) < {krl + kr2So{x)}Sr{x) + kr4, 2 < r < r^,; 
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and 

Vrix) < kr^Sp^r){x) + K^, I < T < T^^.. (2-7) 

are satisfied, where 1 < p(r) < rmix for 1 < r < rmlx- 

The quantities rmL and rmL usually need to be reasonably large, if Assump- 
tion 14.21 below is to be satisfied. 

Now, for Xn as in the introduction, we let denote the time of its 
n-th jump, with = 0, and set := lim^-j-oo ^ , possibly infinite. For 
< t < t^^, we define 

Sl^'Xt) := Sr{X^{t))- Uf\t) := Urix^m V}''\t) := V;(x^(t)), 

(2.8) 

once again with XAr(t) := ^X^it), and also 

r(^)(C) := inf{t < t^^: ^^H^) > NC}, r > 0, (2.9) 

where the infimum of the empty set is taken to be oo. Our first result shows 
that = oo a.s., and limits the expectations of Slf^\t) and S[^\t) for any 
fixed t. 

In what follows, we shall write J^s^^ = a{Xi\f{u),0 < u < s), so that 
{Ts^^ : s > 0) is the natural filtration of the process X^- 

Lemma 2.2 Under Assumptions \2. 1\ t^^ = oo a.s. Furthermore, for any 
t > 0, 

Eisf^t)} < (5f)(0) + Arfco4t)e'=-*; 
E{Si'^)(t)} < {Si''\0) + Nk^4t)e''''. 

Proof. Introducing the formal generator An associated with (11. ip . 

A^fiX) := Nj2MN''X){fiX + J)-fiX)}, X e X+, (2.10) 

Jej 

we note that NUi{x) = AnSi{Nx). Hence, if we define m/^-* by 

Ml''\t) := Si''\t) - Si''\0) -N f U\''\u) du, t > 0, (2.11) 
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for < / < Tmix, it is immediate from f l2.3p . (12 .Sp and (12 .Op that the process 
{M!''\tATi''\C)), t > 0) is a zero mean J-'*^^-' -martingale for each C > 0. 
In particular, considering M['^\t A t[^\C)), it follows in view of (12.61) that 



< Si^^(0)+ r(A;nE{5f)(nArf ^(C))} + iVfci4)rfn. 



Using Gronwall's inequality, we deduce that 

E{Sl''\tA4''\C))} < {S[^\0) + Nh4t)e''''\ (2.12) 
uniformly in C > 0, and hence that 

snp S^{Xn{s)) > NC] < C-^(^i(a;^(0)) + fci4t)e'="* (2.13) 

■0<s<t -I 

also. Hence suPq<^<( S'i(Xjv(s)) < oo a.s. for any t, Mmc^oo Ti^\C) = oo 
a.s., and, from (12. 3 p and (11.30 . it thus follows that = oo a.s. The bound 
on E{5'|^''(t)} is now immediate, and that on E{S'g^''(t)} follows by applying 
the same Gronwall argument to M^^^ {t A (C) ) . I 

The next lemma shows that, if any T > is fixed and C is chosen large 
enough, then, with high probability, sl'^\t) < C holds for all < t < T. 

Lemma 2.3 Assume that Assumptions \2. 1\ are satisfied, and that Sq^\o) < 
NCo and S[^\o) < NCi. Then, for any C > 2(Co + koiT)e^^^^ , we have 

P[{t-o^'^^(C)<T}] < (CiVl)iroo/(iVC2), 
where Kqq depends on T and the parameters of the model. 
Proof. It is immediate from (12.111) and (12.61) that 



< ^^"^^(0)+ [ {koiSif'\u) + Nko^)du+ sup M^^''\u). (2.14) 

Jo 0<u<t 



Hence, from Gronwall's inequality, if Sq^\o) < NCq, then 



S!f\t) < {N{Co + ko^T)+ sup M^">{u)\e 

0<u<t 



Now, considering the quadratic variation of Mq^\ we have 

r /•tAr(^'(C') 

Y,\{Mi''\t^r[''\c'))Y-N 
for any C > 0, from which it follows, much as above, that 
E({Mr(tArr(C^'))}') < E S^N J\i''\u A tI''\C')) du 

< [\kosES[''\u A A''\C')) + Nko,} du. 
Jo 

Using f l2.12p . we thus find that 

E({Mo(^)(tArfHc'))}') < ^iV(Ci + fci4T)(e'^"*-l) + iVfco5t, (2.17) 

uniformly for all C Doob's maximal inequality applied to MQ^\tAT[^\C')) 
now allows us to deduce that, for any C',a > 0, 



koit 



(2.15) 



^0^^^^)^^ ^ = (2.16) 



sup mI,''\uAtI''>{C')) > aN 

■0<u<T 



< 



No? \ku 
say, so that, letting C" — )■ oo 



^ ^^(Ci + A:i4r){e^-"^-l} + fco5T 



02 



sup M^Q^\u) > aN 

■0<u<T 



< 



02 



No? 



also. Taking a = |Ce '^"^^ and putting the result into (12.151) . the lemma 
follows. I 

In the next theorem, we control the 'higher i/- moments' Sr'^\t) of Xjv(t). 
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Theorem 2.4 Assume that Assumptions \2. 1\ are satisfied, and that Si (0) < 
NCi and 5'^JJ)(0) < NC[. Then, for 2 < r < r^AL and for any C > 0, we 
have 

E{Si''\tA4^\C))} < {Si^\0) + Nkr4t)e^''^'^^'^^'^\ 0<t<T. (2.18) 

Furthermore, if for 1 < r < r^L, Si^\o) < NCr and 5jJ^5(0) < NC'^, 
then, for any 7 > 1, 

P[ sup 5f A rl^\C)) > N^C':^] < Krol-^N-\ (2.19) 

0<t<T 

where 

a;^ := (a + Kj" + V(c;. v i))e^^r.+cK2)T 

and Kro depends on C, T and the parameters of the model. 

Proof. Recalling (12. lip , use the argument leading to (I2.12p with the martin- 
gales M^.^\t A t[^\C') a t^^\C)), for any C" > 0, to deduce that 

ESf)(tArr(C")Ar(^)(C)) 

< ^f)(0) + y (^{krl + Ckr2}E[si''\uATi''\C')A4''\C))j+Nkr4)du, 

for 1 < r < TmL, since Sq^\u) < C when u < Tq^\C): define ki2 = 0. 
Gronwall's inequality now implies that 

ESl^\t A 4''\C') A 4''\C)) < {Sl^\0) + Nkr4t)e^''^'^^''^'^\ (2.20) 

for 1 < r < Tmax, and (12.181) follows by Fatou's lemma, on letting C 00. 
Now, also from (12. lip and (12. 6p . we have, for t > and each r < rmL, 

5f)(tAr(^)(C)) 

= 5f )(0) + N / f/W(M) du + M^^\t A 4''\C)) 

Jo 

< Si^\0) + j (^{krl+CK2}Si^\uA4''\C)) + Nkr4^ du 

+ sup Mf)(«Ar(^)(C)). 

0<u<t 
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Hence, from Gronwall's inequality, for all t > and r < rmL, 
Si''\tA4''\C)) < |iV(a + M)+ sup Mf)(MAr('^)(C))|e('=-+^'=-)*. 

0<M<t J 

(2.21) 

Now, as in fl2.16p . we have 

/■iAr(^'(C')Ar^^'(C) 



(2.22) 

from which it follows, using (12. 7p . that, for 1 < r < rmax, 
Ef{Mf)(tArr(C')Arr)(C))}2' 



, .tAr(^){C')Ar(^){C)) 

< f\KsESll]{u A rr(C") A rf )(C)) + iVfc.s} rf« 



this last by fl2.20p . since p(r) < rmL for 1 < r < rmlx- Using Doob's 
inequality, it follows that, for any a > 0, 



P 



sup 

Lo<M<T 



r2 



A^a2 

Taking a = 7a/ (C^ V 1) and putting the result into fl2.2ip gives fl2.19p . with 
Kr^ = {ClKri + Kr2)/{Clyi). ■ 

Note also that suPq<«2. (t) < oo a.s. for all < r < rmL, in view of 
Lemma [2.31 and Theorem 12.41 

In what follows, we shall particularly need to control quantities of the 
form Xljej '^-^(^^('^))'^("^' 0' where x^v := N~^X]^ and 

d{JX) := (2-23) 
i>o 
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for ( & TZ chosen such that ({j) > 1 grows fast enough with j: see fl4.12p . 
Defining 



rW(a,C) := mils: 



J2c^A^N{s))diJX)>ay (2.24) 
Jej J 



infinite if there is no such s, we show in the following corollary that, under 
suitable assumptions, , () is rarely less than T. 

Corollary 2.5 Assume that Assumptions \2. 1\ hold, and that ( is such that 

J2aj{N-^X)d{J,C) < {hN-^Sr{X) + k2}'' (2.25) 
JeJ 

for some 1 < r := r{Q < rmax and some b = b{() > 1. For this value 
of r, assume that Si^\0) < NCr and Sl^J~^{0) < NC^ for some constants 

Cr and C;. Assume further that Sif^\o) < NCo, S[^\o) < NCi for some 
constants Cq, Ci, and define C := 2(Co + kQiT)e^'-^^'^ . Then 

P[r(^)(«,C)<T] < N-\K,olf + K,,{C,yi)C-'}, 

for any a > {^2 + kiC'^j^Y , where 7„ := (a^/^ - / {kiC'^j.} , Kro and C'Jj, 
are as in Theorem \2.4l and Koq is as in Lemma \2.3[ 



Proof. In view of fl2.25p . it is enough to bound the probability 

P[ sup Si''\t) > N{a^/' - k2)/ki]. 

0<t<T 

However, Lemma 12.31 and Theorem 12.41 together bound this probability by 

where is as defined above, as long as a^/'' — k2 > kiC"rp. ■ 

If f l2.25p is satisfied, ^j^j oij{xN{s))d{J, Q is a.s. bounded on < s < T, 

because Sr^\s) is. The corollary shows that the sum is then bounded by 
{^2 + kiC'^.rpY, except on an event of probability of order 0{N~^). Usually, 
one can choose 6 = 1. 
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3 Semigroup properties 

We make the following initial assumptions about the matrix A: first, that 
Aij > for all i ^ j > 0; ^ A^i < oo for all i > 0, (3.1) 

and then that, for some \i G Rj^ such that /i(m) > 1 for each m > 0, and 
for some w > 0, 

< w^/i- (3.2) 

We then use yU to define the /i-norm 

net ^= E'"^"^)!^"! ■■= G 7^: net < oo}. (3.3) 

m>0 

Note that there may be many possible choices for [i. In what follows, it is 
important that F be a Lipschitz operator with respect to the /x-norm, and 
this has to be borne in mind when choosing \i. 
Setting 

where S is the Kronecker delta, we note that Qij > for i ^ j, and that 

using (13.21) for the inequality, so that Qa < 0. Hence Q can be augmented to 
a conservative Q-matrix, in the sense of Markov jump processes, by adding a 
coffin state d, and setting Qig := — ^j>o Qij > 0. Let P(-) denote the semi- 
group of Markov transition matrices corresponding to the minimal process 
associated with Q; then, in particular, 

Q = P'(0) and P'{t) = QP{t) for alH > (3.5) 

(Renter 1957, Theorem 3). Set 

R[jt) := e'"*/i(z)P.,(t)//x(j). (3.6) 
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Theorem 3.1 Let A satisfy Assumptions fl3.ip and fl3.2p . Then, with the 
above definitions, R is a strongly continuous semigroup on TZ^, and 



'^IJ,{i)Rij{t) < /i(j)e"'* for all j and t. 



(3.7) 



i>0 



Furthermore, the sums Ylj>0'^iji^)^jk = {R{t)A)ik are well defined for all 
i, k, and 



A = R'{0) and R\t) = R{t)A forallt>0. 
Proof. We note first that, for x G 71^, 



(3.8) 



\\R{t)x\\, < = e'"*$^$^^(j)^..(t)|: 

«>0 j>0 i>0 j>0 

< e-*5^/i(j)l^,l = e-'\\x\\„ (3.9) 



i>o 



since P(t) is substochastic on Z+; hence R: 71^ — )■ 71^. To show strong 
continuity, we take x G 71^, and consider 



\\R{t)x-x\\^ = ^/i(i) 



i>0 



j>0 



E 

i>0 



I Xi 



< (e-* - 1) 5^5^/.(j)P,.(t)a;, + 5^/x(j)P,.(t)x, + 5^/i(0a;,(l - P..(t)) 

i>0 j>0 i>0 

< - l)J2fiU)x, + 2 5^^(^)x,(l - P,,(t)), 



i>0 



j>0 



i>0 



from which it follows that limf_j,o \\R{t)x — x\\^ = 0, by dominated conver- 
gence, since limt_^o -Pm (^) = 1 for each i>0. 

The inequality (13. 7p follows from the definition of R and the fact that P 
is substochastic on Z+. Then 



fi{i 



/^(i) 



/i(fc) /i(j) 
[(gp(t))., + ^p.,(t)]e-* 



14 



with (QP(t))ij = 'Yl,k>oQikPkj{t) well defined because P{t) is sub-stochastic 
and Q is conservative. Using f l3.5p . this gives 

and this establishes (13.81) . I 



4 Main approximation 

Let Xtv, > 1, be a sequence of pure jump Markov processes as in Section [H 
with A and F defined as in (11. 4p and (II. 5p . and suppose that F: IZ^ — j- 7^^, 
with IZ^ as defined in (13. Sp . for some such that Assumption (13. 2 p holds. 
Suppose also that F is locally Lipschitz in the /x-norm: for any 2; > 0, 

sup ||F(a;)-F(y)|U/||x-y||^ < K{f^,F;z) < 00. (4.1) 
^T^y: I|2:|Im'IIs/IIm<^ 

Then, for x(0) G 71^ and as in (13. 6p . the integral equation 

x(t) = R{t)x{0)+ [ R{t-s)F{x{s))ds. (4.2) 
Jo 

has a unique continuous solution x in 7^^ on some non-empty time interval 
[0,tmax), such that, if tmax < oo, then ||a;(t)||^ — )■ 00 as t — )■ tmax (Pazy 1983, 
Theorem 1.4, Chapter 6). Thus, if A were the generator of R, the function x 
would be a mild solution of the deterministic equations (II. 4p . We now wish 
to show that the process xn '■= N~'^Xn is close to x. To do so, we need a 
corresponding representation for X^. 

To find such a representation, let W{t), t > 0, be a pure jump path on 
that has only finitely many jumps up to time T. Then we can write 

W{t) = W{0) + J2 ^^(^i)' < t < T, (4.3) 

j:aj<t 

where AW{s) := W{s) — W{s—) and aj, j > 1, denote the times when W has 
its jumps. Now let A satisfy (13. ip and (13. 2p . and let R{-) be the associated 
semigroup, as defined in (13. 6p . Define the path W*{t), < t < T, from the 
equation 

W*it) := Rit)WiO) + Zr..,<tRit-^j)^j-IoRit-s)AWis)ds, 

(4.4) 
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where := AW{aj). Note that the latter integral makes sense, because 
each of the sums Ylj>o-^ij{'t)^jk is well defined, from Theorem 13.11 and 
because only finitely many of the coordinates of W are non-zero. 

Lemma 4.1 W*{t) = W{t) for allO<t<T. 

Proof. Fix any t, and suppose that W*{s) = W{s) for all s < t. This is 
clearly the case for t = 0. Let a{t) > t denote the time of the first jump 
of W after t. Then, for any < h < a{t) — t, using the semigroup property 
for R and f O|l . 

W*{t + h) -W*{t) 

= {R{h) - I)R{t)W{0) + J2 (^(^) - ^)^(^ - (4-5) 

j:(7j<t 

ft i-t+h 

- / {R{h)- I)R{t- s)AW{s)ds- I R{t + h - s)AW{t) ds, 
Jo Jt 

where, in the last integral, we use the fact that there are no jumps of W 
between t and t + h. Thus we have 

W*{t + h) -W*{t) 
= {R{h) - I) \ R{t)W{0) + J] - aj)Aj - [ R{t ~ s)AW{s) ds 

/t+h 
R{t + h- s)AW{t) ds 

/t+h 
R{t + h- s)AW{t)ds. (4.6) 



But now, for x G Af. 



t+h 



/t+ii 
R{t + h- s)Ax ds = {R{h) - I)x, 

from ([31]), so that W*{t + h) = W*{t) for all t + h < a{t), implying that 
W*{s) = W{s) for all s < a(t). On the other hand, from (14. 4p . we have 
W*la{t)) - W*{a{t)-) = AW{a{t)), so that W*{s) = W{s) for all s < a{t). 
Thus we can prove equality over the interval [0,cri], and then successively 
over the intervals [aj,aj^i], until [0,T] is covered. I 
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Now suppose that W arises as a realization of X^r. Then X^- has transi- 
tion rates such that 

MN{t) := V AX^ia^)- [ AXN{s)ds- [ NF{xN{s))ds (4.7) 

is a zero mean local martingale. In view of Lemma 14. we can use (14. 4p to 
write 

XN{t) = R{t)XN{0) + MN{t) + N [ R{t- s)F{xN{s))ds, (4.8) 

Jo 

where 

Mjv(t) := R{t- aj)AXNi(rj) 

j:crj<t 

- [ R{t-s)AXN{s)ds- [ R{t-s)NF{xN{s))ds. (4.9) 
Jo Jo 

Thus, comparing (14. 8 p and (14. 2p . we expect xjy and x to be close, for 
< t < T < tmax, provided that we can show that sup^^^ ||mAr(t)||/i is small, 
where m^ij:) := N~^Misr{t). Indeed, if a;Ar(0) and a;(0) are close, then 

\\xNit) - X{t)\\f, 

< ||i?(t)(x^(0)-x(0))t 

+ / ||i?(t-s)[F(x^(s)) -F(x(s))]||^rfs+ ||mjv(t)|U 
Jo 

< e^'\MO)-xm^ 

+ [ e^'-'-'^K{fi,F;2ET)\\xN{s)-x{s)\\^ds + \\mNm^, (4.10) 



by (13.90 ■ with the stage apparently set for Gronwall's inequality, assuming 
that ||xAr(0) — x(0)||^ and supQ<«y ||mAr(t)||^ are small enough that then 
lkAf(i^)IU ^ 2Ht for < t < T, where St := supo<t<T ||2:(^)IU- 



Bounding supo<t<T ||^iv(^)|lAJ is, however, not so easy. Since Mjy is not 
itself a martingale, we cannot directly apply martingale inequalities to control 
its fluctuations. However, since 

Mjv(t) = / R{t~ s)dMN{s), (4.11) 
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we can hope to use control over the local martingale M/v instead. For this 
and the subsequent argument, we introduce some further assumptions. 

Assumption 4.2 

1. There exists r = < rmax such that supj>o{/i(j)/^'r(j)} < c>o. 

2. There exists C, & TZ with C,{j) > 1 for all j such that f l2.25p is satisfied 

(2) 

for some b = b{Q > 1 and r = r{() such that 1 < r{() < rmax, arid that 



The requirement that ( satisfies (I4.12p as well as satisfying fl2.25p for some 
r < rmax implies in practice that it must be possible to take rmix and rmL 
to be quite large in Assumption 12. 1( see the examples in Section O 

Note that part 1 of Assumption 14.21 implies that \imj^^{^{j) / Ur^j)} = 
for some r = < r^ + 1. We define 



where p{-) is as in Assumptions 12. ll We can now prove the following lemma, 
which enables us to control the paths of M^r by using fiuctuation bounds for 
the martingale M^r. 

Lemma 4.3 Under Assumption \4-S^ 




(4.12) 



p(C,/^) 



max{r(C),p(r(C)),r^}, 



(4.13) 



M^it) 



MMit) + [ R{t- s)AMn{s) ds. 



Jo 



Proof. From (13. 8p . we have 



/ + 




Substituting this into (14. lip , we obtain 
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^ 1^ R{v)Al[o,t^s]{v)dv^dMN{s) 



MN{t) 



MN{t) + R{v)Al[o,t-s]iv)dv\dXMis) 



wo 



R{v)Al[o^t-s]{v) dv > Fo{xn{s)) ds. 



It remains to change the order of integration in the double integrals, for 
which we use Fubini's theorem. 

In the first, the outer integral is almost surely a finite sum, and at each 
jump time tf^ we have dX^ltf'^) G JT". Hence it is enough that, for each i, 
m and t, "^j^q Rij{t)Ajm is absolutely summable, which follows from Theo- 
rem [Hill Thus we have 

^ 1^ R{v)Alio^t_4v)dv'^dXN{s) = ^ R{v)A{XN{t~v)-XN{0)}dv. 

(4.14) 

For the second, the k-th component of R{v) AFq{x n (s)) is just 

J2 Rkjiv) ^ A,, Yl J'MMs))- (4.15) 
j>o i>o Jej 

Now, from (13. 7p . we have < Rkj{v) < /i(j)e"'^//i(/c), and 

5^/i(j)|A,,| < iJi{l){2\Au\+w), (4.16) 
i>o 

because A^/i < w^i. Hence, putting absolute values in the summands in fl4.15p 
yields at most 

' Jej z>o 

Now, in view of fl4.12p and since C(j) > 1 for all j, there is a constant K < oo 
such that fi{l){2\Aii\ + w) < K({1). Furthermore, ( satisfies (12.250 . so that, 
by Corollary 12.51 (^j{^n{s)) J2i>o I-^'ICIO is a.s. uniformly bounded in 

< s < T. Hence we can apply Fubini's theorem, obtaining 

R{v)Al[o^t-s]iv)dv\ Fo{xNis))ds = / Riv)A\ Fo{xNis))ds 
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\^J0 



and combining this with fl4.14p proves the lemma. I 

We now introduce the exponential martingales that we use to bound the 
fluctuations of Mat. For 6 G -R^+ bounded and x G 7?.^, 



ZNfi{t) := e^^^'^W exp |- gNe{xN{s-)) rfsj , t> 0, 
is a non-negative finite variation local martingale, where 

gNe[ 



,(0 := $:iVa,(0(e^-^^^^-l). 
For t > 0, we have 



\ogZNfi{t) = 6'^XN{t)- / gNe{xN{s-)) ds 

Jo 

= e^niNit) - ifN,0{^N{s-),s)ds, (4.17) 



where 

meiO ■■= E ^"-^(0 (e"""'"' - 1 - N-'e^j) ^ (4-18) 
Jej 

and mN{t) '■= N~^MN{t). Note also that we can write 

'PnAO = N [\l-r)D\^{^,re)[e,e]dr, (4.19) 

where 

JeJ 

and D^f AT denotes the matrix of second derivatives with respect to the second 
argument: 

I^V(e,^')[Ci,C2] := iV-2E«^(0^''"^''^''^i^^^^C2 (4.20) 

JeJ 

for any (1,(2 e 

Now choose any B := {Bk, > 0) G 7^, and define f^\B) by 



T, 



fi) := inf <^ t > 0: V aj(xjv(t-)) > 5^ 



Our exponential bound is as follows. 
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Lemma 4.4 For any k > 0, 



sup 



0<i<TAf^^'(B) 

for all < 6 < Bi^K^^T , where := J^e'^*, and is as in fll.2p . 

Proof. Take 6 = e'-'^^/?, for /3 to be chosen later. We shall argue by stopping 
the local martingale Z^^g at time a^^\k,S), where 

a'-^'^k^d) := TAfl^\B)Ami{t:m%{t)>5}. 

Note that e^''^^-' < so long as < iV, so that 



\mUt)\>5 



< 2exp{-5'^N/2BkK^T). 



Thus, from (14.191) . we have 



and hence, on the event that a^^\k, 5) = mf{t: m%{t) > 5} < (T Af^^^(S)), 
we have 

ZMA^{k,5)) > exp{(35-^N~'Bk(3^K,T}. 

But since ZN^g{0) = 1, it now follows from the optional stopping theorem 
and Fatou's lemma that 

1 > E{Z^,,(a(^)(fc,5))} 



> P 



sup m%{t) > 6 exp{f36 - \N'^Bkf3^K,T}. 



^0<t<TAf(^'(B) 

We can choose f3 = 6N/BkK^T, as long as 5/BkK^,T < 1, obtaining 



sup m%{t)>5\ < exp{-6'^N/2BkK,T). 



^0<t<TAfl^\B) 

Repeating with 

a^^'^k^d) := T Afjf \B) Ami{t: -m%{t)>6}, 
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and choosing /3 = SN/B^K^T, gives the lemma. I 

The preceding lemma gives a bound for each individual component of M^. 
We need first to translate this into a statement for all components simulta- 
neously. For ( as in Assumption 14.21 we start by writing 

:= maxA:-^#{m:C(m)<A:}; Zf^) := ^M^^^ . (4.21) 

zi^^ is clearly finite, because of Assumption 14.21 and the same is true for zi^^ 
also, since Z of Assumption 14.21 is at least #{m: C{m) < k}/\/k, for each k. 
Then, using the definition f l2.24p of r^^)(a,C), note that, for every fc. 



_(iV) 



for any t < '(a, Q and any h E TZ, and that, for any /C C Z 



+ 5 



'S^ 'S^ ^ ^+\\h^u\ ^ ^ ^ aj{xN{t))h{k)d{JX) 

< • (4.23) 

mmk(.K:{C{k)/h{k)) 

From K22^ with h{k) = 1 for all k, if we choose B := (a/Cik), k > 0), then 
T^^\a, () < fl^\B) for all k. For this choice of B, we can take 

2 .2. ^ 4ai^,TlogAr ABkK^TlogN 

- ^.(«) - iVC(A:) = ^'-'^^ 

in Lemma [4.41 for A; G n^la), where 

fi:^(a) := {A;: ({k) < \aK,TN/logN} = {k: Bt > 4\ogN/K,TN}, 

(4.25) 

since then Sk{a) < Bj^K^T . Note that then, from f l4.12p . 

^ ^^{k)5k{a) < 2ZsJ aK,TN-^\ogN , (4.26) 

with Z as defined in Assumption 14.21 and that 

Ma)\ < laZ^^^K^TN/ log N. (4.27) 
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Lemma 4.5 // Assumptions \4.2\ are satisfied, taking Sk{a) and KAr(a) as 
defined in (I4.24p and fl4.25p . and for any 7] we have 



1. P 



2. P 



U \ sup \mN{t)\ > 6k{a) 



< 



aZi'^K,T_ 
2NlogN' 



J2 X%{t) = for aUO<t<T A r^^^ (a, C) 



> 1 



41ogA^ 



3- sup < y r]{k)\F''{xN{t))\} < . ,,nM n^^■ 



Proof. For part 1, use Lemma 14.41 together with fl4.24p and (14.271) to give 
the bound. For part 2, the total rate of jumps into coordinates with indices 
k ^ KAr(a) is 

a 



^ aj{xN{t)) < 



minfc^«,^(„) C{k)' 



ift < r(^)(a, C), using g;23]) with /C = (/«^(a))^ which, combined with (jOl]) , 
proves the claim. For the final part, if t < T^^\a, (), 

and the inequality follows once more from (I4.23p . ■ 



Let B^^\a) and B\^'{a) denote the events 



>(2), 



k^Kpf(a) 



X%{t) = for all < t < T A T^^\a, Q 



Bf{a) 



Pi \ sup \mN{t)\ < 6k{a) 



(4.28) 



and set B]\i{a) := -B^^(a) fl B)^'{a). Then, by Lemma l475l we deduce that 



?(2) 



P\R (nY] < ^^l!^ , 41ogiV 
- 2iVlogiV 
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(4.29) 



of order 0(A^~^logA^) for each fixed a. Thus we have all the components 
of Mtv simultaneously controlled, except on a set of small probability. We 
now translate this into the desired assertion about the fluctuations of mAr. 

Lemma 4.6 // Assumptions are satisfied, then, on the event Bn^o), 



r- logN 

sup i|m7v(t)IU < y^^iEv 

where the constant depends on T and the parameters of the process. 

Proof. From Lemma [4.31 it follows that 



sup \\mN{t)t (4-30) 

0<i<TAT('V) (a,C) 

< sup ||"^7v(i)IU+ sup / \\R{t - s)Am]sr{s)\\ij_ds. 

0<i<TAT(^)(a,C) 0<i<TAT(^)(a,C) Jo 

For the first term, on -BAr(a) and for < t < T A T^'^\a, (), we have 

The first sum is bounded using fl4.26p by 2Z\JaK^T N^^/'^y/log N, the sec- 
ond, from Lemma [4.51 and ( I4.25p . by 



For the second term in (14.301) . from (13.71) and (14.161) . we note that 
\\R{t ~ s)AmN{s)\\^ < 5^/i(A;)^i?H(t-s)5^|Air||m:^(s)| 

fc>0 l>0 r>0 



Ur\\m''j^[s) 

l>0 r>0 
r>0 
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On Bjv(a) and for < s < TAT^^\a,C), from Km . the sum for r E Kjv(a) 
is bounded using 

J2 Kr){2\A„\+w}\m^^{s)\ 
< Yl f^ir){2\Ar.r\+w}5r{a) 

rGKjv(a) 

'4afs:*TlogA^ 



< l^{r){2\Arr\+w} 



iVC(r) 



< (2V«;)Zv/4air,Tyi^. 



The remaining sum is then bounded by Lemma 14.51 on the set Bn{cl) and 
for < s < T A T^^\a, (), giving at most 

J2 ^i{r){2\Arr\ + w}\m'^{s)\ 



< y\ Kr){2\Arr\+w} [ \F''{xN{t))\dt 



< 



r^K]v(a) 

(2 V w)saJ^, 
{ak)/fi{k){\Akk\ + l}) 



Integrating, it follows that 



sup / \\R{t — s)AmNis)\\^ds 

0<t<TATWia,C) Jo 

I 



lTa\ \okN 



< (2T V l)e"'^ { ^/AaK.TZ + ) J2J, J— ^ 



and the lemma follows. 



This has now established the control on suPq<«2. ||^«Ar(t)||/^ that we need, 
in order to translate (14.101) into a proof of the main theorem. 
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Theorem 4.7 Suppose that (O}, (O}, <^^, and (HH) are all satis- 

fied, and that Assumptions \2. 1\ and \4.^ hold. Recalling the definition ( I4.13P 
of p{(,jj), for ( as given in Assumption \4.^ suppose that S^^^^^^{0) < NC^ 
for some < oo. 

Let X denote the solution to (14. 2 p with initial condition x(0) satisfying 
Spii, I,) ix{0)) < oo. Then tmax = oo. 

Fix any T, and define := supg<«y ||a;(t)||^. // ||xAr(0) — x(0)||^ < 
^'Etc'^^^''*^'^ , where fc* := e'^^K{ix, F; 2Et), then there exist constants c\, c^ 
depending on C^, T and the parameters of the process, such that for all N 
large enough 



sup \\xNit)-x{t)\\,> \e-''\\xNiO)~xm, + ci\M^]e'''' 

0<t<T \ V iV ' 

< (4.31) 

Proof. As 5^2m)(0) < NC,, it follows also that si^\o) < NC, for aU 
< r < p{(, ji). Fix any T < tmax, take C := 2(C* + kQiT)e^°^^ , and observe 
that, for r < p{C,,ld) A rmix, and such that p{r) < p{(,jj), we can take 

Ct < CrT := {2(a V 1) + A;,4T}e('^'-i+^'^'-2)^, (4.32) 

in Theorem 12.41 since we can take to bound Cr and C'^.. In particular, 
r = r{() as defined in Assumption 14.21 satisfies both the conditions on r 
for f02|l to hold. Then, taking a := {fca + kiCrior}^'^^'^ in Corollary EH it 
follows that for some constant C3 > 0, on the event B^la), 

P[r(^n«,C)<r] < CsN-\ 

Then, from (I4.29p . for some constant C4, P[-BAr(a)^] < C4A^~^logA^. Here, 
the constants C3, C4 depend on C^,, T and the parameters of the process. 

We now use Lemma 14.61 to bound the martingale term in (I4.10p . It fol- 
lows that, on the event B^ia) fl {r^^)(a,C) > T} and on the event that 
||a;Ar(s) — a;(s)||^ < Et for all < s < t. 



\xr,{t)-x{t)\\^ < |e"'^||a;^(0)-x(0)||^ + v^i^gj]^^°^^ 



+ k^ I \\xNis) - x{s)\\f,ds, 
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where fc* := e^'^K{^, F;2'E'r). Then from Gronwall's inequahty, on the 
event 5jv(a) n {r(^)(a,C) > T}, 



(4.33) 

for all < t < T, provided that 

This is true for all sufficiently large, if ||xjv(0) - x(0)||^. < |STe-("'+^*)^, 
which we have assumed. We have thus proved (14.311) . since, as shown above, 
P(5iv(a)'^U{r(^)(a,C) > T}") = 0{N-HogN). 

We now use this to show that in fact tmax = oo. For x(0) as above, we 
can take x]v(0) := N-^[Nx^{Q)\ < x^(0), so that ^i^|^)(0) < NC, for : = 

Sp{c,^)ix{0)) < oo. Then, by (|4.13p, limj_,oo{/u(j)/'^p(C,M)0')} = 0' ^o it fol- 
lows easily using bounded convergence that ||xjv(0) — x(0) ||^ — )■ as — )■ oo. 
Hence, for any T < tmax, it follows from (14.311) that ||a;Ar(t) — a;(t)||^ -^d 
as N ^ oo, for t < T, with uniform bounds over the interval, where '— )-_d' 
denotes convergence in distribution. Also, by Assumption 14.21 there is a con- 
stant C5 such that ||a;Ar(t)||^ < c^N^^ Sr'i^\t) for each t, where < rmL and 
< p(C)/^)- Hence, using Lemma [2731 and Theorem 12.41 supg<«2r ||a;Ar(i)||/i 
remains bounded in probability as A^ — )■ 00. Hence it is impossible that 
||x(t)||^ — > 00 as T — )■ tinax < CO, implying that in fact tmax = 00 for such x(0). 




Remark. The dependence on the initial conditions is considerably compli- 
cated by the way the constant C appears in the exponent, for instance in the 
expression for C^t in the proof of Theorem 14.71 However, if kr2 in Assump- 
tions HIT] can be chosen to be zero, as for instance in the examples below, the 
dependence simplifies correspondingly. 

There are biologically plausible models in which the restriction to J' > — 1 
is irksome. In populations in which members of a given type I can fight one 
another, a natural possibility is to have a transition J = —2e^^^ at a rate 
proportional to X\X^ — 1), which translates to aj = a[l^^ = 7a;' (x' — A^~^), 
a function depending on A^. Replacing this with aj = 7(x')^ removes the 
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A^-dependence, but yields a process that can jump to negative values of X . 
For this reason, it is useful to be able to allow the transition rates aj to 
depend on A^. 

Since the arguments in this paper are not limiting arguments for — cx3, 
it does not require many changes to derive the corresponding results. Quan- 
tities such as A, F, Ur{x) and Vr{x) now depend on A^; however, Theorem 14 .71 
continues to hold with constants Ci and C2 that do not depend on A^, provided 
that /i, w, u, the kim from Assumption 12.11 and ( from Assumption 14.21 can 
be chosen to be independent of A^, and that the quantities zi from fl4.2ip 
can be bounded uniformly in A^. On the other hand, the solution x = x'^^^ 
of fl4.2p that acts as approximation to x^v in Theorem 14. 71 now itself depends 
on A^, through R = R^^^ and F = F^^\ If A (and hence R) can be taken 
to be independent of A^, and limAr^oo H-F*-^-* — = for some fixed n~ 
Lipschitz function F, a Gronwall argument can be used to derive a bound 
for the difference between x^^^ and the (fixed) solution x to equation (14. 2 p 
with A^-independent R and F. If A has to depend on A^, the situation is 
more delicate. 

5 Examples 

We begin with some general remarks, to show that the assumptions are sat- 
isfied in many practical contexts. We then discuss two particular examples, 
those of Kretzschmar (1993) and of Arrigoni (2003), that fitted poorly or 
not at all into the general setting of Barbour & Luczak (2008), though the 
other systems referred to in the introduction could also be treated similarly. 
In both of our chosen examples, the index j represents a number of individ- 
uals — parasites in a host in the first, animals in a patch in the second — 
and we shall for now use the former terminology for the preliminary, general 
discussion. 

Transitions that can typically be envisaged are: births of a few parasites, 
which may occur either in the same host, or in another, if infection is being 
represented; births and immigration of hosts, with or without parasites; mi- 
gration of parasites between hosts; deaths of parasites; deaths of hosts; and 
treatment of hosts, leading to the deaths of many of the host's parasites. For 
births of parasites, there is a transition X ^ X + J, where J takes the form 

Jl = 1; Jm = -1; Jj = 0, jj^l,m, (5.1) 
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indicating that one m-host has become an /-host. For births of parasites 
within a host, a transition rate of the form bi^rn'^Xm could be envisaged, 
with I > m, the interpretation being that there are Xm hosts with parasite 
burden m, each of which gives birth to s offspring at rate bg, for some small 
values of s. For infection of an m-host, a possible transition rate would be 
of the form 

since an m-host comes into contact with j-hosts at a rate proportional to 
their density in the host population, and pjr represents the probability of a 
j-host transferring r parasites to the infected host during the contact. The 
probability distributions pj. can be expected to be stochastically increasing 
in j. Deaths of parasites also give rise to transitions of the form f lS.ip . 
but now with I < m, the simplest form of rate being just dmXm for / = 
m — 1, though d = dm could also be chosen to increase with parasite burden. 
Treatment of a host would lead to values of / much smaller than m, and 
a rate of the form nXm for the transition with / = would represent fully 
successful treatment of randomly chosen individuals. Births and deaths of 
hosts and immigration all lead to transitions of the form 

Jl = ±1; Jj = 0, I- (5.2) 

For deaths, Ji = —1, and a typical rate would be d'Xi. For births, Ji = 1, and 
a possible rate would be 'Yl,j>Q^j^'ji (with / = only, if new-born individuals 
are free of parasites). For immigration, constant rates A/ could be supposed. 
Finally, for migration of individual parasites between hosts, transitions are 
of the form 

Jl = Jm = -1; Ji+i = 1; Jm-i = 1; Jj = 0, j ^ l,m,l + i,m-i, 

(5.3) 

a possible rate being 'ymXmN ^Xi. 

For all the above transitions, we can take J* = 2 in (11. 2p . and (II. 3p is 
satisfied in biologically sensible models. (13. ip and (13. 2 p depend on the way in 
which the matrix A can be defined, which is more model specific; in practice, 
(13. ip is very simple to check. The choice of /i in (13. 2 p is influenced by the 
need to have (14. ip satisfied. For Assumptions 12. a possible choice of z/ is to 
take z/(j) = (j + 1) for each j > 0, with Si{X) then representing the num- 
ber of hosts plus the number of parasites. Satisfying (12. 5p is then easy for 
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transitions only involving the movement of a single parasite, but in general 
requires assumptions as to the existence of the r-th moments of the distri- 
butions of the numbers of parasites introduced at birth, immigration and 
infection events. For f l2.6p . in which transitions involving a net reduction 
in the total number of parasites and hosts can be disregarded, the parasite 
birth events are those in which the rates typically have a factor mXm for 
transitions with = —1, with m in principle unbounded. However, at such 
events, an m-individual changes to an m + s individual, with the number s 
of offspring of the parasite being typically small, so that the value of J^Ur 
associated with this rate has magnitude rrf~^] the product mX^m^"^, when 
summed over m, then yields a contribution of magnitude S'r(X), which is al- 
lowable in f l2.6p . Similar considerations show that the terms N~^SQ{X)Sr{X) 
accommodate the migration rates suggested above. Finally, in order to have 
Assumptions 14.21 satisfied, it is in practice necessary that Assumptions 12.11 
are satisfied for large values of r, thereby imposing restrictions on the dis- 
tributions of the numbers of parasites introduced at birth, immigration and 
infection events, as above. 



5.1 Kretzschmar's model 

Kretzschmar (1993) introduced a model 
transitions from state X follows: 

J = e*^*"^^ — e^*^ at rate 

J = — e*^*) at rate 

J = e^^^ at rate 

J = e(*+i) - e(^) at rate 



of a parasitic infection, in which the 

Nifix\ i > 1; 

N{K + ia)x\ i>0; 

NXx'<^{x), i > 0, 



where x := A^~-^X, ip{x) := ||x||ii{c+ ||a;||i}^"'" with c > 0, and ||a;||ii := 
Sj>i il^Pj here, < 6' < 1, and 6^ denotes its i-th power (our 6 corresponds 
to the constant ^ in [7]). Both (11. 2p and (II. 3p are obviously satisfied. For 
Assumptions (13. ip . (13. 2p and (14. ip . we note that equation corresponding 
to ([L5]) has 

A^^ = -{k + i{a + fi)}; Ali_^ = ifi and AJ^ = 0, i > 2; 
An = -{k + a + fi}; A'^q = fi + l39; 
Aqo = -k + (3, i>l, 
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with all other elements of the matrix equal to zero, and 

F'{x) = X{x'-^ -x'Mx), z>l; = -Xx^^{x). 

Hence Assumption (13.11) is immediate, and Assumption (13. 2 p holds for = 
{j + 1)*, for any s > 0, with w = {(3 — For the choice = j + 1, F 
maps elements of TZ^ to 7^^, and is also locally Lipschitz in the /x-norm, with 
K{fi,F;E) = c-^XE{2c + E). 

For Assumptions 12.11 choose = fi; then (12.51) is a finite sum for each 
r > 0. Turning to (12.61) . it is immediate that Uo{x) < (3So{x). Then, for 
r > 1, 

J2MN-'X)X^{{z + 2y-{z + ir} < A|i|||5^rX'(^ + 2)^-1 

j>0 0{ ) ..^g 

since, by Jensen's inequality, Si{X)Sr-i{X) < So{X)Sr{X). Hence we can 
take kr2 = kri = and kri = /3 + r2''"~^X in (12. 6p . for any r > 1, so that 
rL'lx = oo. Finally, for (EZD, 

Vo{x) < {k + l3)So{x) + aSi{x), 

so that ko3 = K + (3 + a and fcos = 0, and 

Vr{x) < r\KS2r{x) + aS2r+l{x) + fiS2r-l{x) + 2^^''-^'> XS2r-l{x)) + f3So{x), 

SO that we can take p{r) = 2r + 1, kr3 = (3 + r'^{K + a + yU + 2'^^^'~^^X}, and 
kr5 = for any r > 1, and so rmax = oo. In Assumptions 14.21 we can clearly 
take = 1 and ({k) = {k + if, giving r(C) = 8, 6(C) = 1 and p(C,/u) = 17. 

5.2 Arrigoni's model 

In the metapopulation model of Arrigoni (2003), the transitions from state X 
are as follows: 

J = e(*-i) - e(*) at rate Nix'idi + 7(1 - p)), i > 2; 

J = e(°) - e(^) at rate Nx\di + 7(1 - p) + k); 

J = e(^+i) - e^^ at rate Nibix\ i > 1; 

J = e(°) - e(^) at rate Nx'k, i > 2; 

J = e^k+i) _ g(fc) ^ g(i-i) _ ^^^g Nix'x^p'^, k>0, i > 1; 
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as before, x := A^~^X. Here, the total number N = J2j>o^j ~ ^oi^) of 
patches remains constant throughout, and the number of animals in any one 
patch changes by at most one at each transition; in the final (migration) 
transition, however, the numbers in two patches change simultaneously. In 
the above transitions, 7,p, k are non- negative, and (di), (bi) are sequences of 
non-negative numbers. 

Once again, both f ll.2p and fll.3p are obviously satisfied. The equation 
corresponding to (11. 4p can now be expressed by taking 

An = -{K + i{bi + di + -f)}] = z(di + 7); = ik, i > I] 

Aoo = -K, 

with all other elements of A equal to zero, and 

F'{x) = p7||a;||ii(x^-^ - x^, i > 1; F\x) = -p-fx°\\x\\u + k, 

where we have used the fact that A^~^ X]j>o"^i ~ Hence Assumption (13. ip 
is again immediate, and Assumption (13. 2p holds for = 1 with w = 0, 
for p(j) = j + 1 with w = maxj(6j — di — •y — k)+ (assuming (6j) and (rfj) 
to be such that this is finite), or indeed for fi{j) = {j + 1)* with any s > 2, 
with appropriate choice of w. With the choice /i(j) = j + 1, F again maps 
elements of TZ^ to TZ^, and is also locally Lipschitz in the p-norm, with 
K{p,F; S) = 3p7S. 

To check Assumptions 12. ![ take u = fi; once again, (12. 5p is a finite sum 
for each r. Then, for (12. 6p . it is immediate that Uq{x) = 0. For any r > 1, 
using arguments from the previous example, 

L i>l i>l k>0 

< r2'''~^{max6i Sr{x) + p'jSi{x)Sr-i{x)} 

i 

< r2'''~^{max6i Sr{x) + p'jSo{x)Sr{x)}, 

i 

SO that, since 5*0(3;) = 1, we can take k^i = r2^~^{maxibi + p7) and kr2 = 
kri = in (12. 6p . and rmL = 00. Finally, for (12.70 . Vq{x) = and, for r > 1, 

Vr{x) 

< <^ 2^^"'^-^^ maxfoj S2r-i{x) + max(i~Mj)S'2r(a;) + 7(1 - p)S2r-i{x) 

+P7(22(^-l)5i(x)52._2(x) + S,{x)S2r-l{x))'] + KS2r{x) , 
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so that we can take p{r) = 2r, and (assuming i ^di to be finite) 
/Cr3 = 1^ + r^{2^''^~^''(max6j + p^y) + max{i^^di) + 7}, 

i i 

and kr5 = for any r > 1, and rmax = 00. In Assumptions 14.21 we can again 
take = 1 and ({k) = {k + 1)'^, giving r{() = 8, 6(C) = 1 and p(C,/^) = 16. 
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